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Abstract 

We study the symmetries of the three heavy-quark system under exchange of the quark fields 
within the effective field theory framework of potential non-relativistic QCD. The symmetries 
constrain the form of the matching coefficients in the effective theory. We then focus on the 
color-singlet sector and determine the so far unknown leading ultrasoft contribution to the static 
potential, which is of order of In fj,, and consequently to the static energy, which is of order of lno s . 
Finally, in the case of an equilateral geometry, we solve the renormalization group equations and 
resum the leading ultrasoft logarithms for the static potential of three quarks in a color singlet, 
octet and decuplet representation. 

PACS numbers: 12.38.-t,12.38.Bx,14.20.-c,14.40.Pq 
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I. INTRODUCTION 



Bound states of a heavy quark Q and an antiquark Q have been the subject of extensive 
theoretical studies since the early days of quantum chromodynamics (QCD). Relatively 
less attention has been paid to bound states of three heavy quarks (QQQ), also referred 
to as triple heavy baryons, as a consequence of their still missing experimental evidence. 
Nevertheless there is an ongoing theoretical activity devoted to their study mostly driven by 



lattice computations |lHl3l|. but also by phenomenologica 



analyses (for a review see [141 ] ) 



and more recently by effective field theory methods [15l4l7l]. The theoretical interest is 
mainly triggered by the geometry of these systems, which allows to address questions that 
are inaccessible with two-body systems. Examples are the minimal energy configuration of 
three quarks in the presence of a confining potential or the origin of a three-body interaction. 
In this paper we will further explore the geometrical properties of the three heavy-quark 
system. 

Systems of heavy quarks are conveniently studied within an effective field theory (EFT) 
framework, a treatment motivated by the observation that these systems are non-relativist ic 
and, therefore, characterized by, at least, three separated and hierarchically ordered energy 
scales: a hard scale of the order of the heavy-quark mass, m, a soft scale of the order of 
the typical relative momenta of the heavy quarks, which are much smaller than m, and 
an ultrasoft (US) scale of the order of the typical binding energy, which is much smaller 
than the relative momenta. 1 We further assume that these scales are much larger than the 
typical hadronic scale Aqcdj i n this way justifying a perturbative treatment for all of them. 
By integrating out modes associated with the different energy scales one goes through a 



sequence o: 
modes 19, 



EFTs 18j|: non-relativistic QCD (NRQCD), obtained from integrating out hard 
20| and potential non-relativistic QCD (pNRQCD), derived from integrating out 

22]. Potential NRQCD provides a formulation 



gluons with soft momenta from NRQCD 



21 



of the non-relativistic system in terms of potentials and US interactions; for this reason 



1 In a thrcc-body system, we may in general expect to have more than one typical relative momentum 
and more than one US energy scale. To keep our discussion simple, we assume all relative momenta to 
be of the same order and so for all US energy scales. In the dynamical case, this is realized when the 
masses of the heavy quarks are of the same order. In the static limit, which will be our main concern in 
the following, this condition is realized by locating the three quarks at distances of the same order. We 
emphasize that this condition may be (also largely) violated in different geometrical configurations. 
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it has proven a convenient framework for calculating US corrections. Although originally 
designed for the study of QQ bound states, i.e. quarkonia, pNRQCD has been subsequently 



15 



16|. 



applied also to baryons with two and three heavy quarks 

In this paper we study the symmetry properties of three heavy-quark systems under 
exchange of the heavy-quark fields and their implications for the form of the pNRQCD 
Lagrangian. We also calculate the US corrections of order In a s to the singlet static 
energy and of order In fi to the singlet static potential of a triple heavy baryon. Whereas 



this has been achieved for the case of QQ systems more than ten years ago [23j, the result 
for QQQ systems will be new. 

The paper is organized as follows. Section [III is devoted to set up pNRQCD for systems 
made of three static quarks. The explicit construction and color structure of the heavy- 
quark composite fields, pNRQCD is conventionally formulated in, is outlined in detail. In 
Sec. [IIH we discuss the symmetry under exchange of the heavy-quark fields and analyze 
its implications for the various matching coefficients, i.e. the potentials, of pNRQCD. In 
Sec. IIVI we determine the correction of order In a s to the singlet static energy. Restricting 
ourselves to an equilateral configuration of the heavy quarks, we finally solve in Sec. [V] 
the renormalization group equations for the singlet, octet and decuplet static potentials at 
leading logarithmic accuracy. We conclude in Sec. [VI] 



II. pNRQCD FOR QQQ 

In this section, we shortly review the basic steps that lead to pNRQCD for systems made 
of three static quarks. Finite mass corrections may be systematically added to the static 



Lagrangian in the form of irrelevant operators, some of which have been considered in [15]. 
The non-relativistic nature of the system ensures that, apart from the kinetic energy, which 
is of the same order as the static potential, 1/m corrections are small. 
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A. NRQCD 



Our starting point is NRQCD in the static limit. In the quark sector the Lagrangian is 



24J and reads 



identical with the heavy-quark effective theory Lagrangian 

£nrqcd = Q j iD°Q + 9 W " ■ (1) 

i 

The heavy-quark fields Q (Q^), which annihilate (create) a heavy quark, are described by 
Pauli spinors, whereas q l are the Dirac spinors that describe light (massless) quarks of flavor /. 
The quantity iD° = id — gA° denotes the time component of the covariant derivative, where 
g is the strong gauge coupling, a s = g 2 /(47r), and A is the time component of the gauge 
field. The Lagrangian (JT|) is insensitive to the flavor assignment of the heavy-quark fields, a 
property known as heavy-quark symmetry. We have omitted the heavy-antiquark sector, as 
it is irrelevant to our scope. 



B. pNRQCD 



For the purpose of studying heavy-quark bound states, it is convenient to employ an EFT 
where the heavy-quark potentials are explicit rather than encoded in dynamical gluons, as 
it is the case in NRQCD. Such an EFT is pNRQCD, which is obtained from NRQCD by 
integrating out gluons whose momenta are soft. The degrees of freedom of pNRQCD are 
heavy-quark fields, light quarks and US gluons. As it is unnecessary to resolve the individual 
heavy quarks, pNRQCD is often formulated in terms of heavy-quark composite fields. The 
matching coefficients of pNRQCD multiplying operators bilinear in the composite fields may 
then be interpreted as the heavy-quark potentials in the corresponding color configurations. 

The derivation of pNRQCD involves identifying the heavy-quark composite fields in 
NRQCD, matching them to pNRQCD, and explicitly ensuring that the resulting pNRQCD 
field content is ultrasoft. We start with the construction of the heavy-quark composite fields. 
This is the point where the specific heavy-quark state that the EFT is meant to describe 
has to be specified. In our case, this is a QQQ state. 
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C. Geometry of a three-quark state 



To characterize the geometry of a QQQ state, we call xi, X2 and X3 the positions of the 
quarks and define the vectors (i — 1, 2, 3) as follows (cf. Fig. [1]), 

ri = xi - x 2 , r 2 = xi - x 3 , r 3 = x 2 - x 3 . (2) 

Note that the three vectors are not independent, for ri +r 3 = r 2 . Moreover, for three quarks 
of equal mass or static, it is useful to define the vectors 




FIG. 1: Triangle formed by three heavy quarks located at the positions xi, x 2 and x 3 . The vector 
A points from the heavy quark at x 3 to the center of mass of the two heavy quarks at xi and x 2 . 

D. Heavy-quark composite fields 

Quarks transform under the fundamental representation, 3, of the (color) gauge group 
SU(3) C . Hence, a generic three (heavy) quark field made of fields located at the same point, 
QiQjQk (hj,k = 1,2,3 denote color indices), transforms as a representation of 3 <8> 3 ® 3. 
The direct product can be decomposed into a sum of irreducible representations of SU(3) C , 
namely 

3®3®3 = 1©8©8©10. (4) 

In general, however, the three quarks are located at different spatial positions Xi, x 2 and 
x 3 . Under an SU(3) C gauge transformation, each heavy-quark field Qi(x, t) transforms as 
<2i(x,t) -»■ Uiii(-x,t)Qii(x,t), where *7(x,t) = exp [i0 a (x, t)T a ], and T a = X a /2 (a = 1,. . .,8) 
denote the eight generators of SU(3) C in the fundamental representation; A a are the Gell- 
Mann matrices. The decomposition requires the fields to be linked to a common point R. 
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For a multi-quark system a natural choice is the system's center of mass. A way to link the 
quark fields to another point is through an equal-time straight Wilson string, 

0(y, x,t) = V exp ji£ jf ds (y - x) • A(x + (y - x)s, t) X , (5) 

where A = A a T a is the color gauge field, and V denotes path ordering of the color matri- 
ces. Due to its transformation property under SU(3) C gauge transformations, 0(y, x, t) — > 
U(y, t)(f>(y, x, t)W(x, t), the Wilson string acts as a gauge transporter, and 0(R, x, t)Q(x, t) 
— > U(R, t)0(R, x, t)Q(x, t) indeed transforms like a quark field located at R. Hence, the 
following three-quark field, 

A^i i/ t(xi,X 2 ,X 3 ,t) = ii /(R,Xi,t)(5^(xi,t)0 J y(R,X2,t)(5jK X 2^)0A : fe'(R,X3,t)Q fc '(x3,t), (6) 

transforms as a 3 <E> 3 <g) 3 representation of the SU(3) C gauge group, and, following Eq. (j4j), 
can be decomposed into a singlet, two octets and a decuplet field with respect to gauge 
transformations in R. 

Since the quark fields do not commute, the order of the quark fields in Eq. ([6]) matters. 
This observation will play a crucial role in Sec. IIIIl For simplicity, we have omitted an explicit 
reference to R in the argument of A4, which includes the time coordinate t and the list of 
position coordinates (xx,X2, X3) of the heavy-quark fields in the order (from left to right) 
of their appearance on the right-hand side of Eq. ([6]). The same convention is used for the 
color indices (i,j,k). 

The composite field Miju may be decomposed into a singlet, S, two octets, A and s , 
and a decuplet, A, according to 

8 

M ijk {xLi, x 2 , xg, t) = S{x.!, x 2 , x 3 , t)S ijk + ° Aa (*i, X 2, x 3 , t)Q$fc 

a=l 

8 10 

+ 3a (xi, x 2, x 3 , t)Q!; k + A 5 (xi, x 2 , x 3 , t)A* jk , (7) 

a=l <5=1 

where S, - fc , O^, O^. and A^ k are orthogonal and normalized color tensors that satisfy the 
relations 

So 1 (~\Aa* (~\Ab jrab r^Sa*f^Sb rah a 8 a 8' cS5' 

—ijk^-ijk — 1 5 ±Ujk ±Ujk ~ ' ±Ujk ±Ujk ~ ' ±±ijk±±ijk — i 

S/-»Aa q r\Sa c a 8 (~\Aa*(~\Sb f~\Aa* a <5 (~\Sa* a 5 r\ /q\ 

—ijk^Lijk — i3.ijk±Lijk — i3.ijk±±ijk — ±Lijk ±Lijk — ±Lijk ±±ijk — ±Lijk ±±ijk ~ u J \°) 
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with a, b G {1, . . . , 8}, and 5, 5' e {1, ... , 10} [15|. If the octet tensors Q^ a k and Of^ have 
the above properties, also the following linear combinations do, 

0$ = (Q# cos W - Og« fe sin w ) , (9) 
Og£ = (0$ sin w + Og, «» w) , (10) 

where u; is an arbitrary angle and ifA, ^Ps denote generic phases. The octet tensors and 
hence form an alternative basis for the 8 © 8 sector. Requiring 

o Aa o Aa k + o Sa o_l a k = o ,Aa a$ + o ,5a og? , (ii) 

the associated octet fields are related to the original ones through the dual relations 

O Mo (x 1 ,x 2 ,x 3 ,0 =e-^[0 Aa (x 1 ,x 2 ,x 3 ,t)cosu;-0 5a (x 1 ,x 2 ,X3,t)sina;] , (12) 
,Sa ( Xl , x 2 , x 3 , t) = e"^ s [0 Aa (x l5 x 2 , x 3 , t) sin w + Sa ( Xl , x 2 , x 3 , t) cos w] . (13) 

To work out the pNRQCD Lagrangian explicitly, we choose a specific (matrix) represen- 



15j, appendix B2. In order to keep this 



tation of the color tensors, namely that given in 
paper self-contained, we reproduce it here. Sticking to this particular choice, the color-octet 
tensors are given by 



Q-tjk — n e ijiKi i (14) 
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and 



Q-ijk - ;r^= { e jkiKi + e iki^ a ji) ■ (15) 

The choice in Eqs. (JHJ) and ffl5|) is such that and O^. are antisymmetric and symmetric 
in the first two color indices, respectively. Consequently, A and s will be referred to as 
the antisymmetric and symmetric octets. Moreover, the color-singlet tensor S ijfc is chosen 
to be totally antisymmetric, 

1 

V6 



Sijjfc - — e ijk , (16) 



whereas the color- decuplet tensor Af- fe is totally symmetric (an alternative decuplet is 



in 



1* 



A 1 - A 4 - A 10 — 1 A 6 
±±111 — ±±.222 ~~ ±±.333 ~~ L ' ±±.{123} 



V6' 



A{112} — A{122} — A{113} — A{223} ~ A{133} ~ A{233} ~ ~f= ■ (l?) 

The symbol {ijk} denotes all permutations of the indices ijk; all components not listed 
explicitly in Eq. (TIT]) are zero. Note that S ijfc and A?- fc are real-valued quantities. 
From Eq. (Q, it follows that the three-quark field 

$jjjfc(xi,X2,x 3 ,t) = Q i (xi,t)Q i (x 2 ,t)(5fe(x3,t) , (18) 

can be written as 

$yfc(xi, x 2 , x 3 , i) = ^«/(xi, R, t)0jj'(x 2 , R, t)0fcfc'(x 3 , R, t)Afiy fc '(xi, x 2 , x 3 , t) , (19) 

where we have used that _1 (y, x, t) = 0'(y, x, t) = 0(x, y, t). Finally, plugging Eq. (J7J) into 
Eq. (fl~9]) we may express the three-quark field $yfc in terms of the composite singlet, octet 
and decuplet fields. The next step will consist in matching these composite fields to the 
corresponding ones in pNRQCD. 



E. Matching and multipole expansion 

We denote with \Q) a generic Fock state containing no heavy quarks, but an arbitrary 
number of US gluons and light quarks: Qi(x, t)\Q) = 0. Therewith we define the three 
heavy-quark Fock state 

\QQQ) = ^J d 3 xx J d 3 x 2 J d 3 x 3 $ ijfc (x 1 ,x 2 ,x 3 ,t)g t fc (x 3 ,t)Qj(x 2 ,t)Qi(x 1 ,t)|fi), (20) 

where M is a normalization factor and the composite field is now interpreted as independent 
of the heavy-quark fields. One can match NRQCD to pNRQCD by equating the expectation 
value of the NRQCD Hamiltonian in the state \QQQ) with the pNRQCD Hamiltonian (see 
18|, |21| for the matching in the QQ case). Thus, the heavy-quark fields in pNRQCD are cast 
into singlet, S, octet, Aa and Sa , and decuplet, A* 5 , fields. The gluons in pNRQCD are 
explicitly rendered US by multipole expanding the gluon fields in the relative coordinates 
rj (i = 1,2,3) with respect to the center of mass coordinate R. The reason is that the 
center of mass coordinate (the "location" of the three heavy-quark system) scales like the 
inverse of the recoiling total momentum of the three quarks, which is of the order of the US 
energy scale, while the relative coordinates of the three quarks (describing the "extension" of 
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the triple heavy baryon) scale like the inverse of the typical relative momenta of the heavy 
quarks, which are of the order of the soft scale. As a result, ultrasoft gluons in pNRQCD 
are invariant under US gauge transformations, i.e. gauge transformations localized in R. 

In the same way as NRQCD can be understood as an expansion of QCD in terms of 
the inverse of the heavy-quark masses, pNRQCD can be understood as an expansion of the 
gluon fields of NRQCD, projected onto the specific (two or three) heavy-quark Fock space, 
in powers of the relative coordinates of the heavy quarks. Quantum corrections of the order 
of the soft scale are encoded in the matching coefficients of pNRQCD in the same way as 
quantum corrections of the order of the heavy-quark masses are encoded in the matching 
coefficients of NRQCD. The matching coefficients of pNRQCD are typically non-analytic 
functions of the relative coordinates. 



F. The pNRQCD Lagrangian 

The pNRQCD Lagrangian is organized as an expansion in 1/m and in the relative co- 
ordinates rj. Up to zeroth order in the 1/m expansion (static limit) and first order in the 



multipole expansion, the pNRQCD Lagrangian for QQQ systems reads [15] 



r - /"(°>°) _i_ /"(o* 1 ) col \ 

''-pNRQCD — ''-pNRQCD "r ''-pNRQCD ■ 

The term ££J^q CD describes at zeroth order in the multipole expansion the propagation 
of light quarks and US gluons as well as the temporal evolution of the static quarks, 
which are cast into singlet, S = 5(xi, x 2 , x 3 , t), octet, A = A (xi, x 2 , x 3 , t) and s = 
O s ('Xi, x 2 , x 3 , t), and decuplet, A = A(x 1; x 2 , x 3 , t), fields (cf. Sec. Ill D| ). 



4£aco = /*d* {* to -V\S + A* [W - V d ] A + O* [iD - K] O* 

+0 5t [iD - V°] s + A] [-VXs\ s + 5t [-VXs\ A } 

+ J]gW -\F; v F a ^ . (22) 
I 

The term £pN^Q CD accounts for the interactions between static quarks and US gluons at 
first order in the multipole expansion, 

4nrqcd = / d 3 P d 3 A { V^ oS f: ^ [5tp • gJSTO*' + O s « p • g&s] 

J a=l 



a=l 



-V, 



O a X-EO a 



- V sx*oA E Te ■ 9^ a Aa + Aa ^X ■ gE a S] 
i 

^ + ^0 A ^X-gE b O' 

,6,c=l 

, [0 Aat p • #E 6 5c + 5at p • gE b O Ac ] 



-V, 



(0,1) 



O s \-EO s 



a,ft,c=l 



a,b,c=l 
8 10 



if abc +3d 



0| 6=1 5=1 
8 10 



a,6=l 5=1 



(23) 



The covariant derivatives, whose time components act on the octet and decuplet fields in 
Eq. f )22|) . are understood to be in the octet and decuplet representations, respectively. They 
are given explicitly in appendix |A] 

A mixing term, — V^ s (0 A ^O s + O s ^O A ), has been included in £^rqcd- Such a term 
was not considered in 15], but was first recognized in [l6|. The mixing potential will play a 
crucial role in the study of the symmetry of pNRQCD under exchange of the heavy-quark 
fields (see Sec. IHIj) and in the calculation of the US corrections to the singlet static energy 
(see Sec. ED- 

For completeness, we list here the leading-order (LO) expressions for the various matching 
coefficients appearing in Eqs. ( |22l) and ( |23l) . At order a s the potentials in Eq. (122]) are given 
by (cf. Q, Q) 



V s (r 1 ,r 2 ,r 3 ) 
V d {r u r 2 ,T 3 ) 
VX(r 1 ,r 2 ,r 3 ) 
^°(ri,r 2 ,r 3 ) 



2 / 1 

3 In 



l 

l r 2| 



1 
l r 3| 



-cc s 



a* 



+ 



+ 



l r l| l r 2| l r 3 

2 1 



11 11 

+ 



3|ri| 12|r 2 | 12 r 3 



1 1 
a s I 77 



5 1 



5 1 



3 rd 12 r 2 12 r 3 



(24) 
(25) 
(26) 
(27) 



10 



V3 



-a* 



1 

l r 2| 



1 

l r 3| 



(28) 



whereas all matching coefficients in Eq. (|23|) are equal to one at LO. The expressions for 
V s up to next-to- next-to-leading order (NNLO), and for V d , V£, Vg and V%s U P to next- 
to-leading order (NLO) can be found in [l6| (the expression for V s up to NNLO is also in 
appendix IB"]) . 



III. SYMMETRY UNDER EXCHANGE OF THE HEAVY-QUARK FIELDS 

As outlined in detail in Sec. Ill D| the heavy-quark fields in the pNRQCD Lagrangian are 
written in terms of composite fields, which are proportional to Qj(x l7 t)Qj(x. 2 , t)Qfc(x 3 , t). 
However, as there is no preferred ordering, and the heavy-quark fields anticommute, different 
orderings of the heavy quarks lead to different composite fields. The orderings are however 
arbitrary and the pNRQCD Lagrangian should be invariant under different orderings of 
the heavy-quark fields. We call this invariance symmetry under exchange of the heavy- 
quark fields or, in short, exchange symmetry. A special case of exchange symmetry is 
the symmetry under permutation of the heavy-quark fields. A different ordering of the 
heavy-quark fields can be realized either (a) by relabeling the heavy-quark coordinates in 
the pNRQCD Lagrangian or (b) by anticommuting the heavy-quark fields in the composite 
fields. Since the two procedures lead to the same Lagrangian, this constrains the form of 
the heavy-quark potentials. In fact, the invariance of the Lagrangian under (a) is trivially 
realized due to the additional integrations over the quark locations x 1; x 2 and x 3 , and only 
(b) results in nontrivial transformations. 

(a) We may relabel the coordinates x, and the relative vectors in the pNRQCD La- 
grangian according to one of the following two possibilities (other relabelings follow from 
these) 



xi O x 2 , x 3 : < 



n ->■ -ri 
r 2 ->■ r 3 , 
r 3 -»■ r 2 

(29) 
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II 



xi -H- x 3 , x 2 : 



ri -»■ -r 3 
r 2 -> -r 2 
r 3 -»■ -ri 



The relabelings affect the pNRQCD potentials and the ordering of the quark fields in the 
composite fields of pNRQCD. 

(b) Because the heavy-quark fields Qi(x) of NRQCD satisfy equal-time anticommutation 
relations, {Qj(x, t), Qj(y,t)} = 0, from Eq. (TT8I) it follows that 



$ ijfc (xi,X2,X 3 ,t) = -$ iifc (x 2 ,Xi,X 3 ,t) 
$ijfe(xi,X 2 ,X3,t) = -$ fcji (x 3 ,X 2 ,Xi,t) 



(30) 
(31) 



These identities hold also for JAijkfei, x 2 , x 3 , t), which is related to (xi, x 3 , x 2 , £) through 
Eq. (USD: 



.W V /,(X|. x 2 . x 3 , t) = -.Mjifc(x 2 , Xi, x 3 , t) , 
Mijk(x-1, x 2 , x 3 , t) = --M fei i(x 3 , X 2 , Xi, t) . 



(32) 
(33) 



In turn, the identities for A^yjfc(xi, x 2 , x 3 , t) enable us to derive corresponding identities for 
the singlet, octet and decuplet fields just by multiplying Eqs. (l3~2~l) and ( |33l) with S ij7 ,, A*- fc , 
or O^f, respectively, and summing over i,j,k: 



and 



5'(xi,x 2 ,x 3 ,t) = 5(x 2 ,xi,x 3 ,t) 
A 5 (xi, x 2 , x 3 , t) = -A <5 (x 2 , xi, x 3 , t) 
Aa (x u x 2 , x 3 , t) = Aa (x 2 , xi, x 3 , t) 
5a (x 1 ,x 2 ,x 3 ,t) = -0 5a (x 2 , Xl) x 3 ,t) 



^(xi, x 2 , x 3 , t) — S*(x 3 , X 2 , Xi, t) 
A 5 (xi, x 2 , x 3 , t) = -A 5 (x 3 , x 2 , xi, t) 

Aa ( Xl , x 2 , x 3 , t) = - cos(f )0 Aa (x 3 , x 2 , Xl> t) + sin(f )0 5a (x 3 , x 2 , Xl , i) 
Sa (x 1; x 2 , x 3 , t) = sin(f )0 Aa (x 3 , x 2 , Xl , t) + cos(f )0 5a (x 3 , x 2 , x 1; t) 



(34) 



(35) 
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At variance with the relabeling (a), anticommuting the heavy-quarks in the composite fields 
only indirectly affects the pNRQCD potentials. Note that the octet transformations in (1341) 
and (|35|) may be interpreted as a special case of the transformations ( fl"2]) and (Tl~3l) . 

By relabeling (a) or by anticommuting the heavy-quark fields (b) we get two versions of 
the pNRQCD Lagrangian that must be the same. This requires the pNRQCD potentials to 
change in a well defined manner under the transformations (i) and (it). In particular, if we 
restrict ourselves to the potentials in Eq. f l22|) . the singlet and decuplet potentials remain 
invariant, whereas the octet potentials transform as 

VZ(-t u t 3 ,t 2 ) =VX(r u r 2 ,v 3 ) 
< V3(-n,r 3 ,r 2 ) =^( ri ,r 2 ,r 3 ) , (36) 
k V As(~ r u r 3 , r 2 ) = -VX s (rx, r 2 , r 3 ) 

and 

VX(-r 3 , -r 2 , -n) = V2(ri, r 2 , r 3 ) cos 2 (f ) + K?(r l5 r 2 , r 3 ) sin 2 (f ) 

-y 4 ° 5 (r 1 ,r 2 ,r 3 )sin(f) 
< V s °(-r 3 , -r 2 , -n) = Vftn, r 2 , r 3 ) sin 2 (f ) + K?(n, r 2 , r 3 ) cos 2 (f ) 

+VX ? (r 1 ,r 2 ,r 3 )sin(f ) 

^ VXs{-t 3 , -r 2 , -n) = \ [VSir,, r 2 , r 3 ) - V^, r 2 , r 3 )] sin(f ) - r 2 , r 3 ) cos(f ) 

(37) 

for transformations of type (i) and (ii) respectively. The above transformation properties 
do not rely on any perturbative expansion; moreover, they are valid also beyond the static 
limit for any order in 1/m. 2 As a simple application of the above formulas, let us consider 
for instance the LO expression of V£g(ri, r 2 , r 3 ) given in Eq. f l28|) . Under fl37j) it transforms 
into 

V^ 5 (-r 3 , -r 2 , -n) = -^a s (±-±V (38) 

4 \l r 2| l r l|/ 

2 Note however that a generalization to finite heavy-quark masses, mi, rri2 and 777.3, would also require 
some adjustment in Eq. ([^1) . as - besides the heavy-quark locations - also the masses would have to be 
exchanged, e.g. in Eq. (J^U) (i), mi ^ m 2 , etc. 
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which is the result as expected from relabeling the coordinates according to Eq. (129]) (ii). 
Let us emphasize again that the inclusion of the octet mixing potential V£ s in Eq. ff22l) is 
essential for reproducing the correct transformation properties of the octet potentials. 

Finally, it is interesting to apply relations f )36|) and f )37|) to the most simple case of an 
equilateral geometry. In such a geometry we have a single length scale r = |i"i| = |r 2 | = |r 3 | 
and a single angle ri • r 2 — • ?3 = ?2 • r 3 = cos (§)- Whenever the potentials are invariant 
under the transformations fl29|) . which is surely the case for two-body interactions but may 
not hold at higher orders, from Eq. (|36|) it follows that V^ s = and from Eq. (|37|) that 

VX(r) = V s °(r) = V°(r) . (39) 



IV. THE QQQ SINGLET STATIC ENERGY AT 0(atlna s 



The potentials of pNRQCD depend in general on a factorization scale [i separating soft 
from US contributions, 3 whereas the singlet static energy E s is an observable and therewith 
independent of \i. As in the QQ case 23j], the QQQ singlet static potential V s is expected 
to become \i dependent at next-to-next-to-next-to leading order (NNNLO), i.e. at order 
Q!g [l||. The difference between the singlet static energy and the singlet static potential is 
encoded in an ultrasoft contribution denoted 5^ s , which starts contributing at order a A s . It 
depends on /x in such a way that E s , given by 

E s (r u r 2 , r 3 ) = V 3 ^, r 2 , r 3 ; //) + ^ s (r l7 r 2 , r 3 ; fi), (40) 

is fi independent. The cancelation of the \x dependence of V s against 5{j S at NNNLO leaves 
in E s a remnant, which is a contribution of order a^\na s . This is the leading perturbative 
contribution to E s that is non- analytic in a s . The most convenient way to calculate the 
In fi term in V s , and the In a s term in E s , is by looking at the leading divergence 
of 5(j S . This requires the one- loop calculation of the color-singlet self energy as opposed to 
the three- loop calculation necessary to extract the term a^ln/i directly from V s . We will 
perform this calculation in the following section. 



3 This dependence, which will be displayed explicitly in the following, has been dropped in Eqs. ([36 
and fl3I]). 
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A. Determination of 5fj S 

We aim at calculating Sfj S up to order a A s . For this purpose we need the singlet and octet 
propagators, and the octet mixing potential at leading order [cf. Eq. (122]) ]. 

= 6(T)e~ (singlet propagator), 

b a = 6(T)e~^sT (symmetric octet propagator), 

" = 9(T)e~^A T S ab (antisymmetric octet propagator), 

(41) 

b x " — b x a = — iV^gSab (octet mixing potential), 



as well as the singlet-to-octet interaction vertices at order in the multipole expansion 
[cf. Eq. ( )23|) . note that the singlet couples differently to the symmetric and antisymmetric 
octets], 




The parameter T in Eq. (HIT) is the propagation time. The wavy lines in Eq. ( )42|) represent 
ultrasoft gluons; note that we have written the vertices with US gluons treating the gluons 
as external fields 

The most noteworthy difference with respect to the calculation of 5(j S in the QQ case is 
that here the singlet couples to two distinct octet fields and that the octet fields mix. For 
this reason the calculation in the baryonic case exhibits some novel features with respect 
to the analogous mesonic case. Since the mixing of the octet fields is an effect of the same 
order as the energies of the octets, it must be accounted for to all orders when computing the 
physical octet-to-octet propagators. The resummation of the octet mixing potential gives 
rise to three different types of resummed octet propagators: 

(1) a resummed octet propagator, G° s , that describes the propagation from a symmetric 
initial state to a symmetric final state: 



15 



(2) a resummed octet propagator, G° A , that describes the propagation from an antisym- 
metric initial state to an antisymmetric final state: 



— E 



n = 



(3) a resummed octet propagator, G AS , that describes the propagation from a symmetric 
initial state to an antisymmetric final state or vice versa: 



The explicit expressions for the resummed octet propagators are most conveniently computed 
in momentum space and read 

~i [G°s(E)} ab = {E _ yo + - v\ hi) - iYXs) 2 ' (43) 

■ \r°(w\\ - t6abi - E ~ v ^ (aa\ 

[ A )lab (E - V s ° + ie){E - V° A + ie) - (V^ 5 ) 2 ' 1 J 

~i [G°As(E)} ab = {E _ yo + te){E ab _ yo + ~ ) _ { yo s) 2 > ( 45 ) 

with e — > + . After performing a Fourier transform from energy E to time T, we obtain 

= *(T) ^ [ (£l " V|) e "^ lT " " ^) ^ ^ ■ 

= 'CO V§) e~^ T - (E 2 - VI) e*&] S ab , (46 ) 



where 
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FIG. 2: Leading-order contributions to 5fj S . As there is no direct coupling between decuplet and 
singlet fields at first order in the multipole expansion, we do not have contributions involving 
decuplet degrees of freedom. 

The US contribution <5{j S is given at LO by the color-singlet self-energy diagrams shown 
in Fig. |2l Because the singlet couples to two distinct octet fields and they mix, we have four 
such diagrams [cf. Eq. f|46]) ]. They give 

-{E 2 - VX)e-^ E2 - ys) ] (p ■ E a (t)p ■ E a (0)> 

-(E 2 - V^e"^ 2 "^] (A ■ E a (t)A ■ E a (0)) 
+2W2 ^7E r ^ ET^E, b~ U{El ~ VS) -e"^- yS) ] (p ■ E a (t)A ■ E"(0)>, (48) 

where (• • ■ ) stands for a vacuum expectation value. In writing the various contributions in 
Eq. f|48l) . we have kept the same order as in Fig. [2j the first two terms correspond to the 
two diagrams shown in the first line of Fig. [21 and the last contribution is the sum of the 
two diagrams in the second line of Fig. [2l which are equal. 

The vacuum expectation value of two chromoelectric fields reads in dimensional regular- 
ization (d = 4 — 2e is the number of dimensions) 

(a • E a (t)b ■ E«(0)> = a • b / ^^M^ + °M , (49) 

where a and b are two generic vectors and t > 0. Performing the integrals in (148]) we obtain 
4a s 1 



6lJS 3 7T Ei - E 2 



- VX) + ^-(Ei - Vs) ~ P ~^V° AS ) {E 1 - V s ) 3 
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1 {E 1 -V s f , 5 
- - 7s - In + - 

e 7r/i z 3 



^(^2 - V° A ) + ^(E 2 - y 5 ) - P -±V° AS ) (£ 2 - v s ^ 3 



(E 2 - V s ) 2 5 



(50) 



where 7^ is the Euler-Mascheroni constant. Equation f )50|) comprises the entire US contri- 
bution up to order a A s . The explicit expressions may be obtained by replacing E\ and E 2 
with the right-hand side of Eq. f|4T|) . and V%, V§ and by the LO expressions given 
Iqs. ( 124]) . (|26|) . (l2Tj) and (128]) respectively. Equation (150]) corrects the expression derived 



111 
in 



151 ] without accounting for the mixing of the octet fields, i.e. assuming V AS = 0. 



B. Invariance of 5fj S under exchange symmetry 

The US correction, Sfj S , calculated in the previous section is expected to be invariant 
under the exchange symmetry discussed in Sec. IIHI To verify this we observe that according 
to Eqs. ([H and ([37} the combinations (VX + V§) and [(V| - Vgf /4 + (V^ 5 ) 2 ] are each 
invariant. This implies that both E\ and E 2 are invariant according to the definition (147]) . 
Also the singlet static potential, V s , is invariant at LO [see Eq. ( 124"]) ]. If we rewrite explicitly 
the expression |p| 2 /4 + |A| 2 /3 in terms of the positions of the heavy quarks with the help of 
Eqs. © and ©, 

— + — = - (x 2 + x 2 + x 2 - Xi ■ X 2 - X! • x 3 - x 2 ■ x 3 ) , (51) 

it is evident that this expression is invariant under the transformations ( l29l) . Finally, we 
have to show that the expression 

Vl^ + V° l -±f + V° AS P -±, (52) 

is also invariant. This is a straightforward, although not manifest, consequence of the trans- 
formations (129]) . (136]) and (]37j) . which completes the proof that <5{j S * s invariant under the 
exchange symmetry. The invariance of Sfj S is directly inherited by the contribution to V s at 
order ln/i and the singlet static energy E s at order lna s . 
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C. The QQQ Singlet Static Potential and Energy 



According to Eq. (|40l) , the divergence and the In /i term in S(j S must cancel against a 
divergence and a term In \x in the singlet static potential V s . Therefore the a* In \x part of 
the potential may be read off from Eq. (150]) . In a minimal subtraction scheme, the singlet 
static potential up to order a\ In /i is then given by 



V s (r 1 ,r 2 ,r 3 ;fi) = V^ NLO (r 1 ,r 2 ,r 3 ) 



■t 1,1,1 
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(53) 

The singlet static potential up to order a^, which we have denoted by V^ NLO , has been 
calculated in Ref. |l6[ and is reproduced in appendix El At order , V^ NLO contains the 
leading three-body potential; also the new term proportional to In fi that we have added 
here is a genuine three-body potential. 

Summing up the singlet static potential (|53|) with the US contribution ( )50|) we obtain the 
singlet static energy up to order a* lna s , which reads 



^(r^r^rg) = V^ NLO (r 1; r 2 , r 3 ) 
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The logarithm of a s signals that an ultraviolet divergence from the US scale has canceled 
against an infrared divergence from the soft scale. 

Finally, it may be useful to express Eqs. ([53]) and (|54l) in a way that makes manifest the 
invariance under exchange symmetry proven in Sec. IIVB1 First, we recall that ri, r 2 and r3 
are not independent (cf. Sec. Ill CP and write 

E s (n, r 2 , r 3 ) = E s (r 2 - r 3 , r 2 , r 3 ) = E s (r 2 , r 3 ), (55) 

then we observe that 

£ s (r 2 ,r 3 ) = £ s (r 3 ,r 2 ). (56) 

Hence an expression of the singlet static energy, which is manifestly invariant under exchange 
symmetry, is 

E-(r u r 2 , r 3 ) = + + E'(-^) (w) 

Similarly one can obtain a manifestly invariant expression of the singlet static potential. 



V. RENORMALIZATION GROUP IMPROVEMENT OF THE SINGLET STATIC 
POTENTIAL IN AN EQUILATERAL GEOMETRY 



The US logarithms that start appearing in the static potential at NNNLO may be re- 
summed to all orders by solving the corresponding renormalization group equations. These 
are a set of equations that describe the scale dependence of the static potentials in the dif- 
ferent color representations. They follow from requiring that the static energies of the QQQ 
system and its gluonic excitations are independent of the renormalization scheme. The po- 
tentials in the different color representations mix under renormalization. This may be easily 
understood by looking at the renormalization group equation for the singlet potential that 
can be derived from /idU s /d/i = — /id<5{j S /d/z and Eq. (15"0]) . 



d T/S 

d/i 



8 Os 
3 71 



v° 



v° 



4 



3 
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o P_± 

y/3 
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V 1 
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IAI 2 



VI? 



+ {v° AS f 
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Yi±n_ v *) 2 +3 <ys-vx)\ 3{VZs) 



(58) 



It shows the explicit dependence of the running of V s on the octet potentials and octet 
mixing potential. 

In the QQ case the renormalization group equations have been solved for the singlet 



static potential at next-to-next-to-leading logarithmic (NNLL]accuracy in 25j and at next 



to-next-to- next-to-leading logarithmic (NNNLL) accuracy in [26|. 4 In the QQQ case similar 
results can be obtained by solving Eq. (|58|) with the corresponding renormalization group 
equations for the octet and decuplet potentials. There is however a difference between 
the QQ and the QQQ case that is worth highlighting. While in a QQ system there is 
just one length, the distance between the heavy quark and antiquark, the generic three- 
body system is characterized by more than one length. For a general three-body geometry, 
therefore, logarithmic corrections in the US scale could be numerically as important as finite 
logarithms involving ratios among the different lengths of the system. The calculation of 
these finite logarithms requires the calculation of the QQQ static Wilson loop. However, 
these logarithms are unimportant if the distances between the heavy quarks are similar. In 
the following, we will therefore restrict ourselves to the simplest case of three static quarks 
located at the corners of an equilateral triangle. In this situation, the three-body system is 
characterized, like the two-body one, by just one fundamental length, which can be identified 
with the length of each side of the triangle: |ri| = |r 2 | = |r 3 | = r. 

In the equilateral limit at least up to NLO, the different octet fields do not mix, moreover, 
as has been shown in Eq. (1391) . the two octet potentials and V% are equal. The US 
contribution for the singlet static energy follows by specializing the general formula fl50l) 
to the equilateral limit. The US contributions for the octet and decuplet static energies 
can be derived along the same lines (cf. also the calculation of the US corrections for the 



QQ octet potential in Ref. [22j). In particular, in the equilateral limit one has to consider 
only the diagrams shown in Fig. [3j since octet-to-octet diagrams with an intermediate octet 
propagator in the loop are scaleless for V§ = V£ = V°, and thus vanish in dimensional 
regularization. Moreover, the US leading-order contribution for the symmetric octet is equal 



4 An NNLL accuracy amounts at resumming a^(a s ln fi) n terms and an NNNLL accuracy amounts at 
resumming a 4 (a s ln^t)™ terms, with n £ Nq. 
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FIG. 3: Leading-order ultrasoft contributions to the singlet, <5{j S , octet, 5^j S , and de- 
cuplet, 5fj S , energies in an equilateral geometry. The triple lines represent the decu- 
plet propagator, 8(T)e~ tVdj Sggr; the decuplet can couple to a symmetric octet, with vertex 
^73 ( e ij'' T ii' T jj'^.'j'k) A-E 6 , or to an antisymmetric octet, with vertex ig (ej^T^Tj^Af, p-E 6 ; 
the other propagators and vertices have been introduced in Eqs. (|4ip and (|42p . 



to the one for the antisymmetric octet; we call it, S^j S . The divergent parts of the diagrams 
shown in Fig. [3] give rise to the following renormalization group equations valid for the 
singlet, octet and decuplet static potentials of three quarks located at the corners of an 
equilateral triangle of side length r: 



du 


= ~«^(V- 


V s ] 


3 + 0{al) 




Sv° 

du 


= ^ [(V - 


-V s 


f + 5{V° - V d f] 




dj. 


= ~ a ,r 2 (V° - 
3n 


yd 


f + O(at) 




d/x 


= a s (3(a s ) 









(59) 



The first equation is just the equilateral limit of Eq. f )58|) . The last equation describes the 
running of the strong coupling constant, where (3(a s ) = — a s /3o/(27r) + 0{a1) is the beta 
function; the first coefficient of the beta function is /?o = H — 2/3n; with ni the number of 
light-quark flavors. By observing that 



V° - V s = -(V° - V d ) + 0{al) 



(60) 



as follows straightforwardly from the results of [16j , the system of equations (159]) can be split 

22 



into two sets of decoupled equations: 
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The two sets of equations can be solved as in [25|] leading to 5 



g o|(l/r) ln « s (l/r) 
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4 /3 r 



a s (/x) 
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In 



(61) 



(62) 



(63) 
(64) 
(65) 



2 (3 r ajji) 

The singlet static potential is known at NNLO, hence Eq. ( )63|) provides the complete ex- 
pression of the singlet static potential at NNLL accuracy in an equilateral geometry. This 
is the most accurate perturbative determination of this quantity. Instead neither the octet 
nor the decuplet potentials are known beyond NLO (see 16|). 



VI. CONCLUSIONS 



In the paper, we have reconsidered the construction of pNRQCD for systems made of 
three heavy quarks with equal masses. We have, in particular, rederived the pNRQCD La- 
grangian in the static limit and put special attention to the symmetry under exchange of 
the heavy-quark fields. Although the symmetry is an obvious property of these systems, its 
consequences for the pNRQCD Lagrangian and in particular for its octet sector have been 
explored here for the first time. Three static quarks may be cast either in a color- singlet, 
two distinct color-octets or a color-decuplet configuration. Whereas the color singlet is com- 
pletely antisymmetric and the color decuplet is completely symmetric in the color-indices, 

5 All coupling constants in V^| NLO (r), VnnloC 7 ") and V^ NLO (r) are evaluated at the scale 1/r. 
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the color-octet transformations depend on the color indices that are exchanged. The fact 
that color-octet fields are specially sensitive to the ordering of the quarks reflects in the fact 
that they mix, in general, under exchange of the heavy-quark fields and dynamically through 
a one-gluon exchange. As a consequence, also the octet potentials and the mixing poten- 
tial transform non trivially under exchange symmetry; we have listed their transformation 
properties in Eqs. (136)) and ( |37|) . 

Thereafter, we have computed the leading ultrasoft contribution to the QQQ singlet static 
energy, 5(j S . Its expression can be found in ( 150]) . Because of the two different octet fields 
and their mixing, the calculation of 5(j S requires the evaluation of four diagrams and the 
resummation of the octet mixing potential for all of them. The calculation is therefore more 
involved than the analogous one of the US contribution in the QQ case. The expression 
for <5{j S in the QQQ case offers also a non-trivial test for the invariance under exchange 
symmetry; this has been performed in Sec. IIVBI A consequence of the calculation of <5{j S 
at leading order is that we can determine the singlet static potential at order at. In /i, see 
Eq. (J53l) . and the singlet static energy at order a$ In a s , see Eq. ( |54l) . These results represent 
the new computational outcome of this work and so far the most accurate determinations 
of the QQQ singlet static potential and energy in perturbative QCD. The new contribution 
computed for the potential is valid for any configuration in space that the three quarks may 
take and it is a three-body interaction. Together with the three-body interaction at two- 



loop order computed in [16J it may provide new insight on the emergence of a long-range 
three-body interaction governed by just one fundamental length that is observed in lattice 
studies (see e.g. {2, 3, 10]). 



In the last part of the paper, we have focused on the special situation where the three 
quarks are located at the corners of an equilateral triangle of side length r. In this limit, 
where the two octet potentials become degenerate, we have solved the renormalization group 
equations for the color singlet, octet and decuplet potentials at NNLL accuracy. The cor- 
responding expressions can be found in Eqs. (|63 l - (j65l) . Hence, for an equilateral geometry, 
the QQQ singlet static potential is now known up to order af(a s \n fj,r) n for all n e Nq. 
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Appendix A: Covariant derivative operators 

In this appendix, we list the explicit matrix representations for the covariant derivative 
operators in the octet and decuplet representations of SU(3) C that appear in Eq. ff22l . The 
SU(3) C covariant derivative is of the general form 

D^ = d, + igA a ^T r a , (Al) 

where a = 1, . . . , 8 and T r a refers to the SU(3) C generators in the representation r. The 



generators in the octet (r = 8) and in the decuplet (r = 10) representation are [15] 
(T«) bc = -tf abc , b,c = l 



)***7 J 

3 
2 



(7? W = o MjkKi'M'jk > S, 5' = 1, . . . , 10, (A2) 



where f abc are the structure constants of SU(3) C . An explicit representation of the decuplet 



tensor Af jk is in (TT7|) . 



Appendix B: The singlet static potential up to order a 



We reproduce here for completeness the expression of the singlet static potential up to 
order computed in [1 61 ] : 



^NNLo( r l' r 2 ? r 3) - -q ^ 
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.3 Vl r i| l r 2| |r 3 
The one-loop and two-loop coefficients a\ and a2 tS depend on the number of light (massless) 
quark flavors, ni, and are given by 
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^ + ^ + l^ + 66((3 )- 484 7 | + 2 04 7£ 
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(B3) 



At two loop, a genuine three-body potential shows up. It is encoded in the function vu 
defined as 
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with R(r 2 ,r 3 ) = xr 2 — yr 3 and M(r 2 ,r 3 ) = |r 2 | yx{l — x) + |r 3 |y^/(I — y). Note that the 
three-body potential in (IB 1 [) is manifestly invariant under the transformations (129)) . 
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